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ABSTRACT 

A spectrum of £ Ori A over the wavelength interval 950— 1150A recorded by 
Interstellar Medium Absorption Profile Spectrograph (IMAPS) on the ORFEUS- 
SPAS I mission shows Lyman and Werner band absorption features from molecular 
hydrogen in rotational levels J = 0, 1, 2, 3 and 5. Most of the molecules are found 



On 
On 

Or 



o 

(N 



> 



in two distinct velocity components. One is at a heliocentric radial velocity of about 
— 1 km s _1 with log A^H^) = 14.5 and a rotational temperature T rot = 950K, while the 
other is at +25 km s" 1 with log A(H 2 ) = 15.9 and T rot = 320K. Some extra H2 exists 



in a much weaker component (log A^H^) = 14.0) between the two main peaks. 



The H2 component at —1 km s 1 exhibits profile shapes that become broader 



0^ ' and show small displacements toward more negative velocities as J increases. These 

changes are inconsistent with a simple interpretation that uv optical pumping in an 
Q\ ■ optically thin, uniform medium creates the H2 in excited rotational levels. Differential 

shielding of the uv radiation at certain velocities does not appear to be a satisfactory 
explanation for the effect. 



Evidence from atomic features at other velocities may offer some insight on the 
origin of this unusual behavior exhibited by the H2 profiles. Absorption features from 
moderately ionized atoms at —94 km s" 1 and more highly ionized species at about 
—36 km s _1 suggest that along the line of sight to £ Ori A there may be a standing 
^ ■ bow shock with an initial compression ratio of 2.6. This shock is probably created 

when a negative-velocity gas flow collides with an obstruction, in this case a neutral 
cloud at km s _1 . If this interpretation is correct, the H2 with the changing profiles 
may represent molecules forming in the postshock gas flow that is undergoing further 
compression as it recombines and cools. We suggest that molecules can form initially 
by associative detachment of H~ in a moving, warm, partly ionized medium behind 
the front. The H2 in this area is most conspicuous in the higher J levels. Later, when 
the gas becomes very cool, neutral, and more compressed as it comes nearly to a halt, 
it is more easily seen in the lowest J levels. In this part of the medium, the principal 
way of producing H2 should be from reactions on the surfaces of dust grains, as one 
expects for quiescent interstellar clouds. 



Subject headings: ISM: molecules — molecular processes — shock waves — stars: 
individual (£ Ori A) 
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1. Introduction 

The first detection of hydrogen molecules in space came from a distinctive pattern of 
absorption features that appeared in a low resolution uv spectrum of £ Per recorded by a 
spectrometer on a sounding rocket (Carruthers 1970). Starting with that pioneering discovery, the 
Lyman and Werner bands of H2 in the spectra of early-type stars have led us down a trail of new 
discoveries about this most abundant molecule in space. Progressively more refined observations by 
the Copernicus satellite have given us a fundamental understanding on this molecule's abundances 
in various diffuse cloud environments (Spitzer et al. 1973; York 1976; Savage et al. 1977), how 
rapidly it is created and destroyed in space (Jura 1974), and the amount of rotational excitation 
that is found in different circumstances (Spitzer & Cochran 1973; Spitzer, Cochran, & Hirshfeld 
1974; Morton & Dinerstein 1976). The observed populations in excited rotational levels have in 
turn led to theoretical interpretations about how this excitation is influenced by such conditions as 
the local gas density, temperature and the flux of uv pumping radiation from nearby stars (Spitzer 
& Zweibel 1974; Jura 1975a, b). Many of the highlights of these investigations have been reviewed 
by Spitzer & Jenkins (1975) and Shull & Beckwith (1982). The Lyman and Werner bands of H2 
can even be used to learn more about the properties of very distant gas systems whose absorption 
lines appear in quasar spectra (Foltz, Chaffee, & Black 1988; Songaila & Cowie 1995), although 
the frequency of finding these H2 features is generally quite low (Levshakov et al. 1992). 

In addition to the general conclusions just mentioned, there were some intriguing details that 
came from the observations of uv absorption lines. The early surveys by the Copernicus satellite 
indicated that toward a number of stars the H2 features became broader as the rotational quantum 
number J increased (Spitzer & Cochran 1973; Spitzer, Cochran, h, Hirshfeld 1974). An initial 
suggestion by Spitzer & Cochran (1973) was that the extra broadening of the higher J levels 
could arise from new molecules that had a large kinetic energy that was liberated as they formed 
and left the grain surfaces. However, a more detailed investigation by Spitzer &; Morton (1976) 
showed that, as a rule, the cases that exhibited the line broadening with increasing J were actually 
composed of two components that had different rotational excitations and a velocity separation 
that was marginally resolved by the instrument. In general, they found that the component with a 
more negative radial velocity was relatively inconspicuous at low J, but due to its higher rotation 
temperature it became more important at higher J and made the composite profile look broader. 

By interpreting the rotational populations from the standpoint of theories on collisional and 
uv pumping, Spitzer & Morton (1976) found a consistent pattern where the components with 
the most negative velocity in each case had extraordinarily large local densities and exposure to 
unusually high uv pumping fluxes. They proposed that these components arose from thin, dense 
sheets of H2-bearing material in the cold, compressed regions that followed shock fronts coming 
toward us. These fronts supposedly came from either the supersonic expansions of the stars' H II 
regions or perhaps from the blast waves caused by supernova explosions in the stellar associations. 

Now, some twenty years after the original investigations with the Copernicus satellite, we 
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have an opportunity to study once again the behavior of the H2 profiles, but this time with a 
wavelength resolution that can cleanly separate the components. We report here the results of 
an investigation of H2 toward £ Ori A, one of the stars studied earlier that showed the intriguing 
behavior with the H2 components discussed above. Once again, the concept of the H2 residing 
in the dense gas behind a shock comes out as a central theme in the interpretation, but our 
description of the configuration given in §5J2 is very different from that offered by Spitzer & 
Morton (1976). 



2. Observations 

The Lyman and Werner band absorptions of H2 in the spectrum of £ Ori A were observed with 
the Interstellar Medium Absorption Profile Spectrograph (IMAPS). IMAPS is an objective-grating 
echelle spectrograph that was developed in the 1980's as a sounding rocket instrument (Jenkins et 
al. 1988, 1989) and was recently reconfigured to fly in orbit. It can record the spectrum of a star 
over the wavelength region 950-1150A at a resolving power of about 200,000.0 This instrument 
flew on the ORFEUS-SPAS carrier launched on 12 September 1993 by the Space Shuttle flight 
STS-51. Jenkins, et al. (1996) have presented a detailed description of IMAPS, how it performed 
during this mission and how the data were reduced. Their article is especially useful for pointing 
out special problems with the data that were mostly overcome in the reduction. It also shows an 
image of a portion of the echelle spectrum of £ Ori A. 

The total exposure time on £ Ori was 2412 s, divided among 63 frames, each of which covered 
74 of the echelle's free spectral range. Spectra were extracted using an optimal extraction routine 
described by Jenkins, et al. (1996), and different measurements of the intensity at any given 
wavelength were combined with weights proportional to their respective inverse squares of the 
errors. Samples of some very restricted parts of the final spectrum are shown in Fig. ffl, where lines 
from J = 0, 1, 3 and 5 may be seen. 



3. Data Reduction 

3.1. Wavelength Scale and Resolution 

Since IMAPS is an objective-grating instrument, there is no way that we can use an internal 
line emission light source to provide a calibration of the wavelength scale. However, as explained in 
Jenkins, et al. (1996), we have an accurate knowledge of how the apparent detector coordinates 
map into real geometrical coordinates on the image plane, and we also know the focal length of 
the cross-disperser grating and the angles of incidence and diffraction for the echelle grating. The 



The observations reported here had a resolution that fell short of this figure, for reasons that are given in 
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Fig. 1. — Sample wavelength intervals covering four H2 lines in the spectrum of ( Ori recorded by 
IMAPS. The intensity scale is in photons per detector pixel. Except for the transition out of J = 5, 
two prominent velocity components are seen in each H2 line, one at a heliocentric radial velocity 
of about —1 km s _1 (labeled with marks that identify the transitions) and another at +25 km s _1 . 
Detector pixels are oversampled by a factor of two in this figure. Dashed lines show the amplitude 
of ten times the expected standard deviation of the points resulting from photon counting statistics. 



- 5 - 



only unknown parameter that we must measure is a zero offset that is driven by the pointing of 
IMAPS relative to the target. We determined this offset by measuring the positions of telluric 
absorption features of O I in excited fine-structure levels. These features are rarely seen in the 
interstellar medium, but there is enough atmospheric oxygen above the orbital altitude of 295 km 
to produce the absorption features in all of our spectra. 

To obtain a wavelength scale that would give heliocentric velocities^] for all of our lines, we 
adjusted the zero offset so that the telluric features appeared at +27.0 km s _1 , a value that 
was appropriate for the viewing direction and time of our observations. The general accuracy 
of our wavelength scale is indicated by the fact that oxygen lines in 4 different multiplets all 
gave velocities within 0.5 km s _1 of the average. Also, for H2 lines out of a given J level that 
had roughly comparable transition strengths, the dispersion of measured velocities was about 
0.5 km s _1 . The measured position of the strongest component (for all J levels) of 24.5 km s _1 
compares favorably with the heliocentric velocity of a strong, but complex absorption feature of 
Na I centered on 24 km s _1 (Welty, Hobbs, & Kulkarni 1994). 

The excited O I lines can also be used to give an indication of the wavelength resolution of 
our observations. We measured equivalent widths of 10.5 and 7.25mA for the O I* and O I** lines 
at 1040.9 and 1041. 7A, respectively. For the applicable densities and temperatures of the Earth's 
upper atmosphere, the occupation of the singly excited level (O I*) should be 3 times that of 
the doubly excited level (O I**), i.e., their relative numbers are governed by just their statistical 
weights g. Making use of this fact allows us to apply a standard curve of growth analysis to derive 
log N(0 I**) = 14.19 and b = 0.99 km s -1 (equivalent to a doppler broadening for T = 950K).f] 
The observed profiles have widths that correspond to b = 3.0 km s _1 , which leads to the conclusion 
that the instrumental spread function is equivalent to a profile with 



(The representative b for the excited O I lines has been elevated to 1.25 km s 1 to account for the 
small broadening caused by saturation). 

The wavelength resolving power that we obtained is lower than what is achievable in principle 
with IMAPS and the pointing stability of the spacecraft. We attribute the degradation to small 
motions of the echelle grating during the exposures, caused by a sticky bearing that relieved 



2 To obtain the LSR velocity in the direction of £ Ori A, one should subtract 17.5 km s _1 from the heliocentric 
velocity. Differential galactic rotation at an assumed distance of 450 pc to £ Ori A should cause undisturbed gases in 
the general vicinity of the star to move at 4.5 km s _1 with respect to the LSR if the galaxy has a constant rotation 
velocity of 220 km s" 1 and Ro = 8.5 kpc (Gunn, Knapp, & Tremaine 1979). Thus, any feature appearing at a 
heliocentric velocity of 22.0 km s _1 should be approximately in the rest frame of gaseous material in the vicinity of 
our target. 

3 These results agree very well with predictions of the MSIS-86 model of the Earth's thermosphere (Hedin 1987) 
for the column density and temperature along a sight line above our orbital altitude and at a moderate zenith angle 



V3.0 2 - 1.25 2 = 2.7 km s 
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(40°). 
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mechanical stresses at random times. The magnitude and character of this effect is discussed in 
detail by Jenkins, et al. (1996). 



3.2. Absorption Line Measurements 

We used the MSLAP analysis program^ to define the continuum level Jo and re-express the 
intensities I(v) in terms of the apparent absorption optical depths function of radial 

velocity, 

T.M = ln( J | J ) . (2) 

For the ideal case where the instrument can resolve the finest details in velocity, T a (v) usually 
gives an accurate depiction of a differential column density per unit velocity through the relation 

N a (v) = 3.768 x 10 14 ^^cm- 2 (km s" 1 )" 1 , (3) 

where / is the transition's /-value and A is expressed in A. However, if there are saturated, 
fine-scale details that are not resolved, the true optical depths t(v) averaged over velocity will 
be underestimated, and one will miscalculate the true column density N(v). One can ascertain 
that this is happening if the application of Eq. ^ for weaker lines indicates the presence of more 
material than from the strong ones (Savage & Sembach 1991; Jenkins 1996). As will be evident 



in §4.1, this appears to happen for the strongest features of H2 in the J = 0, 1 and 2 levels of 



rotational excitation. 

For J levels through 3, we were able to draw together the results for many different 
absorption lines, each going to different rotational and vibrational levels in the upper electronic 
states, 2p<7.B 1 £^ and 2pir C 1 T1 U . In so doing, it was important to keep track of the errors in the 
measured I(v) and combine redundant information at each velocity in a manner that lowered the 
error in the final result. To achieve this goal, we evaluated for every individual velocity point the 
X 2 from a summation over the separate transitions, 

The expected errors in intensity e(J/Jo) represented a combination of several sources of error: 
(1) the noise in the individual measurements of /, (2) an error in the placement of the continuum 
Iq, and (3) an error in the adopted value of zero spectral intensity (which is a finite value of real 
intensity extracted from the echelle order). The errors in / (item 1) were measured from the 
dispersion of residual intensities on either side of the adopted continuum at points well removed in 



4 MSLAP is a third-generation program developed for NASA. MSLAP is copyrighted by Charles L. Joseph and 
Edward B. Jenkins. 
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velocity from the absorption feature. This error generally becomes larger at progressively shorter 
wavelengths, because the sensitivity of IMAPS decreases. (Variations of sensitivity also result 
from being away from the center of the echelle blaze function.) In every case, the noise errors 
were assumed to be the same magnitude at low I/Iq at the centers of lines because statistical 
fluctuations in the background illumination are important. (Generally, the background was about 
as large as Io, so the noise amplitude would decrease only by a factor of y2.) In a number of 
cases, the computed S/N was higher than 50 (see Tables HJ— @). Because there might be some 
residual systematic errors that we have not recognized, we felt that it was unwarranted to assume 
that these cases had the full reliability as indicated by the calculation of S/N, when compared 
with other measurements at lower S/N. To account for this, we uniformly adopted an estimate for 
the relative noise level consistent with the value 



The error in Io (item 2 in the above paragraph) represents the uncertainty of the continuum 
level that arises from a pure vertical translation that would be permitted by the noise in the many 
intensity measurements that define Io. It does not include errors in the adopted curvature of the 
continuum [see a discussion of this issue in the appendix of Sembach & Savage (1992)]. For most 
cases, the curvature was almost nonexistent. The error in the adopted background level (item 3) 
was judged from the dispersion of residual intensities of saturated atomic lines elsewhere in the 
spectrum. At every velocity point, the worst combinations of the systematic errors (i.e., both the 
adopted continuum and background levels are simultaneously too high or, alternatively, too low) 
were combined in quadrature with the random intensity errors (item 1), as modified in Eq. |5], to 
arrive at the net e{I /Iq). 

Tables p]— [| show the transitions for the four lowest rotational levels of H2 covered in our 
spectrum of £ Ori. Laboratory wavelengths are taken from the calculated values of Abgrall, et al. 
(1993a) for the Lyman band system and Abgrall, et al. (1993b) for the Werner bands. Transition 
/-values are from Abgrall & Roueff (1989). The listed values of S/N are those computed as 
described above, but without the modification from Eq. 

All of the lines for J = 4 were too weak to measure. Only one line from J = 5 was strong 
enough to be useful (the Werner 0-0 Q(5) line at 1017.831A with log(/A) = 1.39), although 
weaker lines showed very noisy profiles that were consistent with this line. 

Many lines (or certain portions thereof) were unsuitable for measurement. These lines and 
the reasons for their rejection are discussed in the endnotes of the tables. Table ^ omits some lines 
that are far too weak to consider. 




(5) 
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Table 1. Lines from J=0 



Ident. a 


A (A) 


Log ( fX) 


S/N 


0-0 R(0) b 


1108.127 


0.275 


50 


1-0 R(0) 


1092.195 


0.802 


77 


2-0 R(0) 


1077.140 


1.111 


46 


3-0 R(0) 


1062.882 


1.282 


45 


4-0 R(0) 


1049.367 


1.383 


30 


5-0 R(0) 


1036.545 


1.447 


39 


6-0 R(0) 


1024.372 


1.473 


36 


7-0 R(0) 


1012.810 


1.483 


81 


8-0 R(0) 


1001.821 


1.432 


32 


9-0 R(0) c 


991.376 


1.411 




10-0 R(0) d 


981.437 


1.314 


23 


11-0 R(0) e 


971.985 


1.289 




12-0 R(0) d 


962.977 


1.098 


14 


13-0 R(0) d 


954.412 


1.126 


20 


W 0-0 R(0) f 


1008.552 


1.647 


31 


W 1-0 R(0) s 


985.631 


1.833 




W 2-0 R(0) h 


964.981 


1.823 





a All transitions are in the 2per B X S+ «— X 1 E+ Lyman band system, unless preceded with 
a "W" which refers to the 2p7rC 1 IL i <~ X 1 ^ Werner bands. 

b Not used in the composite profile, because components 1 and 2 were too weak compared 
with the noise. For component 3, this was the weakest line and had the least susceptibility 
to errors from saturated substructures. This line was used to define the preferred value for 
Ntatal with Method A (see § |4.1.l|) . 

c Not considered because this line had interference from the W 1—0 P(3) line. 

d Not included in the composite profile because the S/N was significantly inferior to those 
of other lines of comparable log (/A). 

°Stellar flux severely attenuated by the Ly-7 feature. 
f Not included; there is serious interference from the W 0—0 R(l) line. 
g Not included; there is serious interference from the W 1—0 R(l) line. 
h Not included; there is serious interference from the W 2—0 R(l) line. 
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Table 2. Lines from J=l 



Ident. a 


A (A) 


Log (/A) 


S/N 


0-0 P(l) b 


1110.062 


-0.191 


46 


1-0 P(l) c 


1094.052 


0.340 


40 


2-0 P(l) 


1078.927 


0.624 


33 


3-0 P(l) 


1064.606 


0.805 


48 


4-0 P(l) 


1051.033 


0.902 


48 


5-0 P(l) 


1038.157 


0.956 


78 


6-0 P(l) d 


1025.934 


0.970 




7-0 P(l) 


1014.325 


0.960 


62 


8-0 P(l) 


1003.294 


0.931 


19 


9-0 P(l)° 


992.808 


0.883 


14 


10-0 P(l)° 


982.834 


0.825 


14 


11-0 P(l)° 


973.344 


0.759 


3 


12-0 P(l)° 


964.310 


0.683 


12 


13-0 P(l)° 


955.707 


0.604 


9 


W 0-0 Q(l) 


1009.770 


1.384 


36 


W 1-0 Q(l) f 


986.796 


1.551 


8 


W 2-0 Q(l) f 


966.093 


1.529 


10 


0-0 R(l) s 


1108.632 


0.086 


39 


1-0 R(l) h 


1092.732 


0.618 


69 


2-0 R(l) 


1077.700 


0.919 


55 


3-0 R(l) 


1063.460 


1.106 


59 


4-0 R(l) 


1049.960 


1.225 


64 


5-0 R(l) 


1037.149 


1.271 


56 


6-0 R(l) 1 


1024.986 


1.312 


11 


7-0 R(l) 


1013.434 


1.307 


48 


8-0 R(l) 


1002.449 


1.256 


16 


9-0 R(l) j 


992.013 


1.252 


19 


10-0 R(l) c 


982.072 


1.138 


14 


11-0 R(l) c 


972.631 


1.134 


5 


12-0 R(l) c 


963.606 


0.829 


12 


13-0 R(l) c 


955.064 


0.971 


9 
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Table 2 — Continued 



Ident. a 


A (A) 


Log (/A) 


S/N 


W 0-0 R(l) k 


1008.498 


1.326 




W 1-0 R(l) 1 


985.642 


1.512 




W 2-0 R(l) m 


965.061 


1.529 





a All transitions are in the 2per B X S+ «— X 1 E+ Lyman band system, unless preceded with 
a "W" which refers to the 2p7rC 1 II„ <- X 1 !^ Werner bands. 

b For component 3, this was the weakest line and had the least susceptibility to errors 
from saturated substructures. This line was used to define the preferred value for iVtotal- 
Component 1 of the 0—0 R(2) is near this feature, but it is not close and strong enough to 
compromise the measurement of iVtotal with Method A (§ [4.1.1[ ). We did not use the line in 
the composite profile however. 

c Not used in the composite profile because of interference from the 1—0 R(2) line. 
This interference did not compromise our use of the line for obtaining a measurement of 
Component 3 using Method B (§ [1.1. 2j ). 

d Stellar flux severely attenuated by the Ly-/3 feature. 

°Not included in the composite profile because the S/N was significantly inferior to those 
of other lines of comparable log (/A). 

f S/N too low to use this line, even though its log(/A) is large. 

g Components 1 and 2 too weak to measure, hence not included in composite profile. For 
Component 3, this line was used in Method B (§4.1.2). 

h Possible interference from 1092.620 and 1092. 990A lines of S I, hence not included in 
composite profile. 

'On a wing of the stellar Ly-/3, hence the S/N is low. Line not used in the composite 
profile. 

JNot included in the composite profile because the error array shows erratic behavior. 
k This line has interference from the W 0—0 R(0) and 8—0 P(3) lines. It was not used. 
'This line has interference from the W 1—0 R(0) line. It was not used. 
m This line has interference from the W 2—0 R(0) line. It was not used. 



- 11 - 



Table 3. Lines from J=2 



lUcIlX . 


A [A) 


T orr ( f \\ 


c; /"vr 

u f IN 


0-0 P(2) b 


1112.495 


-0.109 


39 


1-0 P(2) c 


1096.438 


0.420 


51 


2-0 P(2) 


1081.266 


0.706 


55 


3-0 P(2) d 


1066.900 


0.879 


65 


4-0 P(2) 


1053.284 


0.982 


35 


5-0 P(2) 


1040.366 


1.017 


38 


6-0 P(2)° 


1028.104 


1.053 


13 


7-0 P(2) f 


1016.458 


1.007 


34 


8-0 P(2) 


1005.390 


0.998 


29 


9-0 P(2)s 


944.871 


0.937 


18 


10-0 P(2) h 


984.862 


0.907 


5 


11-0 P(2) h 


975.344 


0.809 


7 


12-0 P(2) h 


966.273 


0.798 


13 


13-0 P(2) h 


957.650 


0.662 


12 


W 0-0 P(2) ; 


1012.169 


0.746 


23 


W 1-0 P(2) h 


989.086 


0.904 


8 


W 2-0 P(2) h 


968.292 


0.843 


14 


W 0-0 Q(2) 


1010.938 


1.385 


29 


W 1-0 Q(2)J 


987.972 


1.551 


7 


W 2-0 Q(2)J 


967.279 


1.530 


11 


0-0 R(2) k 


1110.119 


0.018 


45 


1-0 R(2) ! 


1094.243 


0.558 


56 


2-0 R(2) c 


1079.226 


0.866 


35 


3-0 R(2) 


1064.994 


1.069 


53 


4-0 R(2) 


1051.498 


1.168 


76 


5-0 R(2) 


1038.689 


1.221 


72 


6-0 R(2) m 


1026.526 


1.267 




7-0 R(2) 


1014.974 


1.285 


52 


8-0 R(2) 


1003.982 


1.232 


40 


9-0 R(2) 


993.547 


1.228 


20 


10-0 R(2) 


983.589 


1.072 


18 


11-0 R(2) h 


974.156 


1.103 


4 


12-0 R(2) n 


965.044 


0.161 




13-0 R(2) h 


956.577 


0.940 


10 


W 0-0 R(2) 


1009.024 


1.208 


32 


W 1-0 R(2) h 


986.241 


1.409 


3 
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Table 3 — Continued 



Ident. a 


A (A) 


Log (/A) 


S/N 


W 2-0 R(2) 


965.791 


1.490 


16 



a All transitions are in the 2pcr B <— X 1 S+ Lyman band system, unless preceded with 
a "W" which refers to the 2 P tt C x U u <- X X S+ Werner bands. 

b Component 1 of this line is too weak to see above the noise, and Component 3 has 
interference from Component 1 of the 0—0 R(3) line. Hence this transition is not useful. 

c In constructing the composite profile, we used only the velocity interval covering 
Component 3 because Components 1 and 2 arc completely buried in the noise. 

d Component 1 feature seems to be absent for some reason that is not understood. Perhaps 
an unidentified feature on the edge of this component makes it unrecognizable. 

e Stellar flux severely attenuated by the Ly-/3 feature. This line was not used because the 
S/N was too low. 

f Component 1 was badly corrupted by an unidentified line. Only the region around 
Component 3 was used. 

g The nearby W 1—0 Q(5) line makes the continuum uncertain. Thus, we did not use the 
9-0 P(2) line. 

h Not included in the composite profile because the S/N was significantly inferior to those 
of other lines of comparable log (/A). 

'This line was not used because it might be corrupted by the presence of the 1012. 502A 
line of S III at -80 km s _1 . 

j S/N too low to use this line, even though its log(/A) is large. 

k We used only the portion covered by Component 3, since Component 1 of this line has 
serious interference from Component 3 of the 0—0 P(l) line. 

'We used only the portion covered by Component 3, since the continuum just to the left 
of Component 1 is compromised by the presence of Component 3 of the 1—0 P(l) line. 

m Stellar flux severely attenuated by the Ly-/3 feature. This line was not used. 
n There is interference from the N I line at 965.04lA. Hence this line was not used. 
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Table 4. Lines from J=3 



Trlonf a 
ICltJIlI. 


A [A) 




O / IN 


0-0 P(3) b 


1115.895 


-0.083 


45 


1-0 P(3) b 


1099.787 


0.439 


31 


2-0 P(3) c 


1084.561 


0.734 




3-0 P(3) 


1070.141 


0.910 


26 


4-0 P(3) 


1056.472 


1.006 


56 


5-0 P(3) 


1043.502 


1.060 


48 


6-0 P(3) 


1031.192 


1.055 


41 


7-0 P(3) 


1019.500 


1.050 


57 


8-0 P(3) d 


1008.383 


1.004 




9-0 P(3) c 


997.824 


0.944 


18 


10-0 P(3)° 


987.768 


0.944 


10 


11-0 P(3)° 


978.217 


0.817 


20 


12-0 P(3)°' f 


969.089 


0.895 


10 


13-0 P(3)° 


960.449 


0.673 


12 


W 0-0 P(3) 


1014.504 


0.920 


54 


W 1-0 P(3) g 


991.378 


1.075 




W 2-0 P(3) c 


970.560 


0.974 


10 


W 0-0 Q(3) 


1012.680 


1.386 


31 


W 1-0 Q(3) h 


989.728 


1.564 




W 2-0 0(3)^ 


969.047 


1.530 


8 


0-0 R(3) k 


1112.582 


-0.024 




1-0 R(3) L 


1096.725 


0.531 




2-0 R(3) 


1081.712 


0.840 


47 


3-0 R(3) 


1067.479 


1.028 


42 


4-0 R(3) 


1053.976 


1.137 


37 


5-0 R(3) 


1041.157 


1.222 


49 


6-0 R(3) 


1028.985 


1.243 


24 


7-0 R(3) 


1017.422 


1.263 


35 


8-0 R(3) 


1006.411 


1.207 


18 


9-0 R(3) 


995.970 


1.229 


33 


10-0 R(3) c 


985.962 


0.908 


4 


11-0 R(3) m 


976.551 


1.104 




12-0 R(3) c 


967.673 


1.347 


10 


13-0 R(3) c 


958.945 


0.931 


10 


W 0-0 R(3) n 


1010.129 


1.151 


40 


W 1-0 R(3)' 


987.445 


1.409 


6 
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Table 4 — Continued 



Ident. a 


A (A) 


Log (/A) 


S/N 


W 2-0 R(3) c 


966.780 


0.883 


16 



a All transitions are in the 2ptr B : E+ <— X 1 S+ Lyman band system, unless preceded with 
a "W" which refers to the 2 P tt C 1 n„ <- X X S+ Werner bands. 

b This line is too weak to show up above the noise. It was not used in constructing the 
composite profile. 

c This line could not be used because it has serious interference from the 1084.562 and 
1084. 580A lines from an excited fine-structure level of N II. 

d Not used since this line has interference from the W 0—0 R(l) line. 

e Not included in the composite profile because the S/N was significantly inferior to those 
of other lines of comparable log (/A). 

f Not used since this line has interference from the W 2— 0Q(3) line. 

s Not used since this line has interference from the 9—0 R(0) line. 

h Line is submerged in a deep stellar line of N III at 989. 8A. Thus, it could not be used. 

'S/N too low to use this line, even though its log (/A) is large. 

JNot used since this line has interference from the 12— 0P(3) line. 

k This line could not be used because it has interference from the 0—0 P(2) line. 

'This line could not be used because it has interference from the 1096. 877A line of Fe II. 

m The left-hand side of Component 1 has interference from Component 3 of the line of O I 
at 976.448A. This line could not be used even for Component 3 because the continuum level 
was uncertain. 

"Line inadvertently omitted. The omission was discovered long after the combined analysis 
had been completed. 
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4. Results 

Figs. Q—^ show gray-scale representations of x 2 — X m in as a function of logJV (t;) and the 
heliocentric radial velocity v. The minimum value Xmin was determined at each velocity, and 
our representation that shows how rapidly x 2 increases on either side of the most probable 
logN a (v) (i.e., the value where Xmin * s achieved) is a valid measure of the relative confidence of 
the result (Lampton, Margon, & Bowyer 1976). Since we are measuring a single parameter, the 
X 2 distribution function with 1 degree of freedom is appropriate, and thus, for example, 95% of 
the time we expect the true intensity to fall within a band where x 2 ~ Xmin < 3.8, i.e., the "±2<t" 
zone. To improve on the range of the display without sacrificing detail for low values of x 2 — Xmin' 
the actual darknesses in the figures and their matching calibration squares on the right are scaled 
to the quantity log(l + X 2 ~ Xm in )- Measurements at velocities separated by more than a single 
detector pixel (equivalent to 1.25 km s" 1 ) should be statistically independent.^ This separation 
is less than the wavelength resolving power however. Thus, reasonable assumptions about the 
required continuity of the profiles for adjacent velocities can, in principle, restrict the range of 
allowable departures from the minimum x 2 even further than the formal confidence limits. 

The profiles that appear in Figs. ||— 1| indicate that there are two prominent peaks in H2 
absorption, with the left-hand one holding molecules with a higher rotational temperature than 
the one on the right. This effect, one that creates dramatic differences in the relative sizes of the 
two peaks with changing J, was noted earlier by Spitzer, et al. (1974). There is also some H2 that 
spans the velocities between these two peaks. For the purposes of making some general statements 
about the H2, we identify the material that falls in the ranges —15 to +5, +5 to +15, and +15 to 
+35 km s _1 as Components 1, 2 and 3, respectively. While some residual absorption seems to 
appear outside the ranges of the 3 components, we are not sure of its reality. Some transitions 
seemed to show convincing extra absorption at these large velocities, while others did not. 

Component 1 shows a clear broadening as the profiles progress from J = to 5. Precise 
determinations of this effect and the accompanying uncertainties in measurement will be presented 



in §[4^. The widths of the profiles for Component 3 also seem to increase with J, but the effect is 
not as dramatic as that shown for Component 1. We are reluctant to present any formal analysis 
of the broadening for Component 3 because we believe the N a (v) profile shapes misrepresent the 



true distributions of molecules with velocity for J = 0, 1 and 2, for reasons given in §4.1. As a 
rough indication of the trend, we state here only that the apparent profile widths are 4.5, 5.8, 5.8 
and 7.7 km s _1 (FWHM) for J = 0, 1, 2 and 3, respectively. These results for this component 
only partly agree with the finding by Spitzer, et al. (1974) that the velocity width of molecules 
in the J = 1 state is higher than those in both J = or J = 2. The latter conclusions were based 
on differences in the b parameters of the curves of growth for the lines. 



5 This statement is not strictly true, since single photoevents that fall near the border of two pixels will contribute 
a signal to each one. However, the width of one pixel is a reasonable gauge for distance between nearly independent 
measurements if one wants to judge the significance of the x 2 ' s - 
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Fig. 2. — A composite of 8 absorption profiles from H2 in the J = rotational level. Transitions 
listed in Table |l| were used, except where noted. Shades of gray, as indicated by the boxes on the 
right, map out the changes in \ 2 ~ Xmin as a function of \ogN a {v) for each value of v. For reasons 
discussed in §4.1 the strong peak on the right-hand side probably understates the true amount of 
H2 that is really present. 
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Fig. 3. — Same as for Fig. ^, except that the applicable transitions, listed in Table are from the 
J = 1 level. Thirteen transitions were combined to make this figure. As with Fig. ||, the rightmost, 
strong peak probably under-represents the true amount of H2. 
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Fig. 4. — Same as for Fig. ^, except that the applicable transitions, listed in Table ||, are from the 
J = 2 level. A total of 19 transitions were used to construct this figure, but because of interference 
problems only 14 of them covered Components 1 and 2. 
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Fig. 5. — Same as for Fig. 0, except that the 15 applicable transitions, listed in Table ||, are from 
the J = 3 level. Unlike the cases for J = 0, 1 or 2, the right-hand peak does not show any disparities 
in the height from one transition to another, indicating that the representation is probably correct. 




Fig. 6 



. — Same as for Fig. ^, except that only one transition, the Werner 0—0 Q(5) line, was used. 
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Table || lists our values for the column densities / N a (v) dv, obtained for profiles that follow 
the valley of minimum \ 2 - Exceptions to this way of measuring A^H^) are discussed in §4.1 below. 
We also list in the table the results that were obtained by Spitzer, et al. (1974) and Spitzer 
& Morton (1976). With only two significant exceptions, our results seem to be in satisfactory 
agreement with these previous determinations. One of the discrepancies is the difference between 
our determination log N( J = 5) = 13.70 for Component 1, compared with the value of 13.32 found 
by Spitzer, et al. (1974). We note that latter was based on lines that had special problems: 
either the lines had discrepant velocities or the components could not be resolved. The second 
discrepancy is between our value of log N(J = 0) = 15.09 (Method A discussed in § |4.1.1|) or 14.79 
(Method B given in §4.1.2) for Component 3 and the value 15.77 found by Spitzer &; Morton 
(1976) from an observation of just the Lyman 4— 0R(0) line. However, this line is very badly 
saturated (the central optical depth must be about 12 with our value of log iV and 6 = 2 km s _1 ), 
and thus it is not suitable, by itself, for measuring a column density. 



4.1. Unresolved Saturated Substructures in Component 3 

For the right-hand peaks (Component 3) in J = 0, 1 and 2, the weakest transitions show more 
H2 than indicated in Figs. ^ to ||, which are based on generally much stronger transitions. This 
behavior reveals the presence of very narrow substructures in Component 3 that are saturated and 
not resolved by the instrument. Jenkins (1996) has shown how one may take any pair of lines (of 
different strength) that show a discrepancy in their values of N a (v), as evaluated from Eqs. [2] and|||, 
and evaluate a correction to r a (v) of the weaker line that compensates for the under-representation 
of the smoothed real optical depths t(v). In effect, this correction is a method of extrapolating 
the two distorted N a (v)'s to a profile that one would expect to see if the line's transition strength 
was so low that the unresolved structures had their maximum (unsmoothed) t(v) <C 1. 

Unfortunately, we found that for each of the three lowest J levels, different pairs of lines 
yielded inconsistent results. In each case, an application of the analysis of the first and second 
weakest lines gave column densities considerably larger than the same procedure applied to the 
second and third weakest lines. We list below a number of conjectures about the possible cause(s) 
for this effect: 

1. The functional forms of the distributions of subcomponent amplitudes and velocity widths 
are so bizarre, and other conditions are exactly right, that the assumptions behind the 
workings of the correction procedure are not valid. As outlined by Jenkins (1986, 1996), 
these distributions would need to be very badly behaved. 

2. We have underestimated the magnitudes of the errors in the determinations of scattered 
light in the spectrum, which then reflect on the true levels of the zero-intensity baselines 
and, consequently, the values of T a {v) near maximum absorption. 
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Table 5. Log Column Densities 3, and Rotational Temperatures 





Component 1 


Component 2 b 


Component 3 


J 


(-15 < v < +5 km s" 1 ) 


(+5 < v < +15 km s -1 ) 


(+15 < v < +35 km s -1 ) 





13.53 (13.46B? 13.48*) 


12.86 (13.23 d ) 


15.09? 14.79 f (15.2lA? 15.77 d ) 


1 


13.96 (14.15BJ 14.20*) 


13.44 (13.85*) 


15.72? 15. 69 f (15.43B c ) 


2 


13.64 (13.64B? 13.68*) 


13.27 (13.36*) 


14.78? 14.66 f (14.74B? 14.87*) 


3 


13.99 (14.05A? 14.08 d ) 


13.55 (13.69 d ) 


14.19 (14.14A? 14.34 d ) 


4 


••• (13.22A? 13.11 d ) 


••• (-^o d ) 


••• (12.95A? 12.85 d ) 


■5 


13.70 (13.32A? 13.45 d ) 


13.13 g (12.48 d ) 


13.21 (<12.79?-oo d ) 


Total 


14.52 


14.01 


15.86? 15. 79 f 


Rot. Temp. h 


950K 


960K 


320K? 340K f 



a Numbcrs in parentheses arc from earlier Copernicus results reported by Spitzer, Cochran & 
Hirshfeld (1974) and Spitzer & Morton (1976) for comparison with our results (and to fill in for J 
= 4). 

b Not really a distinct component, but rather material that seems to bridge the gap between 
Components 1 and 3. 

c From Spitzer, Cochran & Hirshfeld (1974), with errors A = 0.04-0.09 and B = 0.10-0.19. 

d From Spitzer & Morton (1976). 

Derived from Method A discussed in § [4.1.l| . 

f Derived from Method B discussed in § |4.1.2| . 

s Not a distinct component (see Fig. |^). The number given is a formal integration over the specified 
velocity range and represents the right-hand wing of the very broad component centered near the 
velocity of Component 1. 

h From the reciprocal of the slope of the best fit to \n[N(J)/g(J)] vs. Ej, excluding J = 4. 
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3. The transition /-values that we have adopted are wrong. The sense of the error would be 
such that the weakest lines are actually somewhat stronger than assumed, relative to the 
/-values of the next two stronger lines. Another alternative is that the second and third 
strongest lines are much closer together in their /-values than those that were adopted. 

While we can not rigorously rule out possibilities (1) and (2) above, we feel that they are 
unlikely to apply. Regarding possibility (3), the /-values are the product of theoretical calculations, 
and to our knowledge only some of the stronger transitions have been verified experimentally (Liu 
et al. 1995). It is interesting to see if there is any observational evidence outside of the results 
reported here that might back up the notion that alternative (3) is the correct explanation. 

We are aware of two potentially useful examples where the weakest members of the Lyman 
series have been seen in the spectra of astronomical sources. One is in a survey of many stars by 
Spitzer, Cochran & Hirshfeld (1974) ,f] and another is an array of H2 absorption features identified 
by Levshakov & Varshalovich (1985) and Foltz, et al. (1988) at z = 2.811 in the spectrum of 
the quasar PKS 0528—250. The quasar absorption lines have subsequently been observed at much 
higher resolution by Songaila & Cowie (1995) using the Keck Telescope. 

In the survey of Spitzer et al. (1974), the only target that showed lines from J=0 that were 
not on or very close to the flat portion of the curve of growth (or had an uncertain measurement 
of the Lyman 0—0 R(0) line) was 30 CMa. The 10mA equivalent width measured for this line 
is above a downward extrapolation of the the trend from the stronger lines. If the line's value 
of log /A were raised by 0.28 in relation with the others, the measured line strength would fall 
on their adopted curve of growth. Unfortunately, we can not apply the same test for the Lyman 

0— P(l) or 0— R(2)lines, the two weakest lines that we could use here for the next higher J 
levels, because these lines were not observed by Spitzer et al. 

The H2 lines that appear in the spectrum PKS 0528—250 are created by a heavy-element gas 
system that is moving at only 2000 to 3000 km s~ x with respect to the quasar (and hence one that 
is not very far away from the quasar). The overall widths of the H2 lines of about 20 km s _1 were 
resolved in the R = 36,000 spectrum of Songaila & Cowie (1995), but the shallow Lyman 0—0, 

1— and 2— 0R(0) features showed a strengthening that was far less than the changes in their 
relative /-values. Songaila & Cowie interpreted this behavior as the result of saturation in the 
lines if they consisted of a clump of 5 unresolved, very narrow features, each with b = 1.5 km s _1 , 
distributed over the observed velocity extent of the absorption. One might question how plausible 
it is to find gas clouds with such a small velocity dispersion that could cover a significant fraction 
the large physical dimension of the continuum-emitting region of the quasar. As an alternative, 
we might accept the notion that the lines do not contain unresolved saturated components, but 
instead, that the real change in the /-values is less than assumed. 



6 There are many papers that report observations of H2 made by the Copernicus satellite. Oddly enough, the 
paper by Spitzer et al. (1974) is the only one that includes measurements of the weakest lines. 
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Finally, we turn to our own observations. In our recording of the Lyman 0—0 R(0) line in our 
spectrum of £ Ori, the amplitude of the T a (v) profile of Component 1 (about 4a above the noise), 
in relation to that of Component 3, is not much different than what may be seen in the next 
stronger line, 1— 0R(0). If significant distortion caused by unresolved, saturated substructures 
were occurring for Component 3 in the latter, the size difference for the two components would be 
diminished, contrary to what we see in the data. If one were to say that the difference in log(/A) 
for the two lines were smaller by 0.4, we would obtain N a (v)'s that were consistent with each 
other. 

We regard the evidence cited above as suggestive, but certainly not conclusive, evidence 
that our problems with the disparity of answers for N a (v) might be caused by incorrect relative 
/-values. Even if this conjecture is correct, we still do not know whether the stronger or weaker 
/-values need to be revised. In view these uncertainties, we chose to derive N(v) for Component 3 
by two different methods, Method A and Method B, outlined in the following two subsections. 
Total column densities / N(v) dv derived each way are listed in Table [|. 



4.1.1. Method A 

Method A invokes the working assumption that the adopted /-value for the weakest line is 
about right, and that there is a problem with the somewhat stronger lines. If this is correct, then 
our only recourse is to derive N a (y) from this one line through the use of Eqs. Q and |^ and assume 
that the correction for unresolved saturated substructures is small. For J = 0, 1 and 2, we used 
the Lyman 0—0 R(0), 0—0 P(l) and 0—0 R(2) lines, respectively. (The weakest line for J = 2, 
0— OP (2) could not be used; see note b of Table ||[) 



4.1.2. Method B 

Here we assume that the published /-value for the weakest line is too small, but that the 
values for the next two stronger lines are correct. We then derive corrections for r a {v) for the 
weaker line using the method of Jenkins (1996). While the errors in this extrapolation method 
can be large, especially after one considers the effects of the systematic deviations discussed earlier 



[items (2) and (3) covered in £3^2], under the present circumstances they are probably not much 
worse than the arbitrariness in the choice of whether Method B is any better than Method A or 
some other way to derive N a (v). Lyman band line pairs used for this method were 1—0 R(0) and 
2-0 R(0) for J = 0, 0-0 R(l) and 1-0 P(l) for J = 1, and 1-0 P(2) and 1-0 R(2) for J = 2. 
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4.2. Profile Changes with J for Component 1 

Figures ^ to [6| show very clearly that the profiles for Component 1 have widths that 
progressively increase as the rotational quantum numbers go from J = to J = 5. Figure [7| shows 
a consolidation of the results from Figs. || to ^: the valleys of \ 2 — Xmm are depicted as lines [now 
in a linear representation for iV a (H2)], and the profiles are stacked vertically to make comparisons 
for different J in Component 1 more clear. In addition to showing the changes in profile widths, 
this figure also shows that there is a small (~ 1 km s _1 ) shift toward negative velocities with 
increasing J up to J = 3, followed by a more substantial shift for J = 5. 

A simple, approximate way to express numerically the information shown in Fig. [?] is to 
assume that most of the H2 at each J level has a one-dimensional distribution of velocity that 
is a Gaussian function characterized by a peak value for N(v), iV max , a central velocity, vq, and 
a dispersion parameter, b. We can then ascertain what combinations of these 3 parameters give 
an acceptable fit to the data as defined, for example, by values \ 2 ~ Xmm < ^-8 that lead to a 
95% confidence limit. We carried out this study with x 2 ' s > 01 the type displayed in Figs. [2] to |(| 
summed over velocity points spaced 1.6 km s _1 apart to assure statistical independence. Table |6] 
summarizes the results of that investigation. The quantities v m m and f max are the velocity limits 
over which the fits were evaluated. The error bounds are defined only by the \ 2 limits and do not 
include systematic errors, such as those that arise from errors in /-values or our overall adopted 
zero-point reference for radial velocities. For given J levels, there are small differences between the 
preferred log(A r max ^/7r6) and the log column densities given in Table [5] caused by real departures 
from the Gaussian approximations (J = 5 shows the largest deviation, 0.08 dex, as one would 
expect from the asymmetrical appearance shown in Fig. [?]). 

To determine the real widths of the profiles, one must subtract in quadrature two sources of 
broadening in the observations. First, there is the instrumental broadening of each line in the 



spectrum that we recorded, as discussed in §3.1. Adding to this effect are the small errors in 



registration of the lines, as they are combined to create the x 2 ~ Xmm plots (Figs. |2| to |6|). From 
the apparent dispersion of line centers at a given J, we estimate the rms registration error to be 
0.5 km s _1 . We estimate that the effective b parameter for these two effects combined should be 
about 2.8 km s _1 , and thus the formula given in note a of Table |6| should be applied to obtain a 
best estimate for the true b of each H2 profile (the results for the lowest J levels will not be very 
accurate, since 6 Q b s is only slightly greater than 2.8 km s^ 1 ). 

The results shown in Fig. |7| and Table || show two distinct trends of the profiles with increasing 
J. First, the most probable values for the widths b increase in a steady progression from J = 
to J = 5. Second, the most probable central velocities vq become steadily more negative with 
increasing J, except for an apparent reversal between J = 2 and J = 3 that is much smaller than 
our errors. It is hard to imagine that systematic errors in the observations could result in these 
trends. The absorption lines for different J levels appear in random locations in the spectral image 
formats, so any changes in the spectral resolution or distortions in our wavelength scale should 
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Fig. 7. — Plots of log N a (v) versus v scaled such that the heights of the peaks for Component 1 are 
nearly identical. To facilitate comparisons of widths and velocity centers across different J levels, 
the two vertical, dotted lines mark the half amplitude points of the J = profile. 
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Table 6. Gaussian Fits to Component 1 





J = 


J = 1 


J = 2 


J = 3 


J = 5 


f max (km s 1 ) 


+4 


+4 


+5 


+5 


+7 


^min 


-5 


-6 


-8 


-8 


-14 


Largest log iV max [cm _2 (km s -1 ) -1 ] 


12.82 


13.12 


12.68 


13.00 


12.66 


Most probable logiV max 


12.77 


13.09 


12.62 


12.96 


12.56 


Smallest logiV max 


12.72 


13.06 


12.58 


12.94 


12.46 


Largest vq (km s'^ 1 ) a 


-0.3 


-0.9 


-1.0 


-1.0 


-1.0 


Most probable 


-0.5 


-1.0 


-1.5 


-1.3 


-2.9 


Smallest 


-0.7 


-1.2 


-2.0 


-1.6 


-4.4 


Largest b (km s _1 ) b 


3.2 


4.2 


7.0 


6.8 


14 


Most probable b 


2.9 


3.9 


6.0 


6.5 


9.4 


Smallest b 


2.6 


3.8 


5.2 


6.0 


7.2 



a Heliocentric radal velocity of the profile's center. 

b Includes instrumental broadening and registration errors (see §4.2). Hence, the 

^2 



real b should equal about yb^ — (2.8 km s 



-28- 



affect all J levels almost equally. 

5. Discussion 

5.1. Preliminary Remarks 

The information given in Table [| shows that the 3 molecular hydrogen velocity components 
toward Q Ori A have populations in different J levels that, to a reasonable approximation, conform 
to a single rotational excitation temperature in each case. This behavior seems to reflect what 
has been observed elsewhere in the diffuse interstellar medium. For instance, in their survey of 28 
lines of sight, Spitzer, et al. (1974) found that for components that had N(J = 0) <; 10 15 cm~ 2 , a 
single excitation temperature gave a satisfactory fit to all of the observable J levels. By contrast, 
one generally finds for much higher column densities that there is bifurcation to two temperatures, 
depending on the J levels [see, e.g., Fig. 2 of Spitzer & Cochran (1973)]. This is a consequence 
of the local density being high enough to insure that collisions dominate over radiative processes 
at low to intermediate J and thus couple the level populations to the local kinetic temperature, 
whereas for higher J the optical pumping can take over and yield a somewhat higher temperature. 
For cases where the total column densities are exceptionally low [7V(H2) ~ 10 13 cm -2 for such 
stars as C Pup, 7 2 Vel and r Sco], the rotation temperatures can be as high as about 1000K. This 
behavior is consistent with what we found for our Components 1 and 2. Our Component 3 has a 
somewhat lower excitation temperature, but one that is in accord with other lines of sight that 
have 7V(H2) ~ 10 15 cm~ 2 in the sample of Spitzer, et al. (1974). 

It is when we go beyond the information conveyed by just the column densities and study 
changes in the profiles for different J that we uncover some unusual behavior. Here, we focus 
on Component 1, where the widths and velocity centroids show clear, progressive changes with 
rotational excitation. While Component 3 also shows some broadening with increasing J, the 
magnitude of the effect is less, and it is harder to quantify because there are probably unresolved, 
saturated structures that distort the N a (v) profiles. The changes in broadening with J are 
inconsistent with a simple picture that, for the most diffuse clouds, the excitation of molecular 
hydrogen is caused by optical pumping out of primarily the J = and 1 levels by uv starlight 
photons in an optically thin medium. 

We might momentarily consider an explanation where the strength of the optical pumping 
could change with velocity, by virtue of some shielding in the cores of some of the strongest 
pumping lines. However, in the simplest case we can envision, one where the light from £ Ori 
dominates in the pumping, the shielding is not strong enough to make this effect work. For 
example, in Component 1 we found log./V(H2) = 13.53 (Table |B|) and a largest possible real value^\ 



7 See note b of Table | 
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of b = 1.55 km s~ for molecules in the J = level. We would need to have a pumping line from 
J = with a characteristic strength log /A = 3.0 to create an absorption profile 1 — I(v)/Iq that is 
saturated enough to have it appear, after a convolution with our instrumental profile, as broad as 
the observed N a (v) for molecules in the J = 2 state.Q In reality, the strongest lines out of J = 
have log fX only slightly greater than 1.8 (see Table [j]). Likewise, the width of the N a (v) profile 
for J = 3 can only be matched with a pumping line out of J = 1 with log /A = 2.0, again a value 
that is much higher than any of the actual lines out of this level (see Table ^|). Thus, if we are to 
hold on to the notion that line shielding could be an important mechanism, we must abandon the 
idea that £ Ori is the source of pumping photons. 

We could, of course, adopt a more imaginative approach and propose that light from another 
star is responsible for the pumping. Then, we could envision that a significant concentration of H2 
just off our line of sight could be shielding (at selective velocities) the radiation for the molecules 
that we can observe. While this could conceivably explain why the profiles for J > 1 look different 
from those of J = or 1, it does not address the problem that the profile for J = 1 disagrees with 
that of J = 0. (The coupling of these two levels by optical pumping is very weak.) As indicated 
by the numbers in Table ^, both the velocity widths and their centroids for these lowest two levels 
differ by more than the measurement errors. 

Another means for achieving a significant amount of rotational excitation is heating due to 
the passage of a shock — one that is slow enough not to destroy the H2 (Aannestad & Field 
1973). Superficially, we might have imagined that Component 1 is a shocked portion of the gas 
that was originally in Component 3, but that is now moving more toward us, relatively speaking. 
However this picture is in conflict with the change in velocity centroids with J, for the gas would 
be expected to speed up as it cools in the postshock zone where radiative cooling occurs. Our 
observations indicate that the cooler (rear) part of this zone that should emphasize the lower J 
levels is actually traveling more slowly. 

From the above argument on the velocity shift, it is clear that if we are to invoke a shock 
as the explanation for the profile changes, we must consider one that is headed in a direction 
away from us. If this is so, we run into the problem that we are unable to see any H2 ahead of 
this shock, i.e, at velocities more negative than Component 1. Thus, instead of creating a picture 
where existing molecules are accelerated and heated by a shock, we must turn to the idea that 
perhaps the molecules are formed for the first time in the dense, compressed postshock zone, out 
of what was originally atomic gas undergoing cooling and recombination. In this case, one would 
look for a shock velocity that is relatively large, so that the compression is sufficient to raise the 
density to a level where molecules can be formed at a fast rate. 



8 This simple proof is a conservative one, since it neglects other processes that tend to make the J = 2 profile as 
narrow as that for J = 0, such as pumping from many other, much weaker lines or the coupling of molecules in the 
J = 2 state with the kinetic motions of the gas through elastic collisions. 
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5.2. Evidence of a Shock that could be Forming H2 

5.2.1. Preshock Gas 

There is some independent evidence from atomic absorption lines that we could be viewing 
a bow shock created by the obstruction of a flow of high velocity gas coming toward us, perhaps 
a stellar wind or a wind-driven shell (Weaver et al. 1977). A reasonable candidate for this 
obstruction is a cloud that is responsible for the low-ionization atomic features that can be seen 
near v = km s _1 . 

In the IMAPS spectrum of £ Ori A, there are some strong transitions of C II (1036. 337A) and 
N II (1083. 990A) that show absorption peaks at —94 km s _1 , plus a smaller amount of material at 
slightly lower velocities (Jenkins 1995). Features from doubly ionized species are also present at 
about the same velocity, i.e., C III, N III, Si III, S III (Cowie, Songaila, & York 1979) and Al III 
(medium resolution GHRS spectrum in the HST archive^). Absorption by strong transitions of 
O I and N I are not seen at —94 km s _1 however. The moderately high state of ionization of this 
rapidly moving gas, a condition similar to that found for high velocity gas in front of 23 Ori by 
Trapero et al. (1996), may result from either photoionization by uv radiation from the Orion 
stars or collisional ionization at a temperature somewhat greater than 10 4 K. 

Figure |8| shows spectra that we recovered from the HST archive^ in the vicinity of the C IV 
doublet (1548.2, 1550.8A). We determined an upper limit logiV(C IV) < 12.2 at v « -90 km s _1 . 
When this result is compared with the determination logA r (C II) = 13.84 (Cowie, Songaila, 
& York 1979) or 13.82 (IMAPS spectrum), we find that T < 20,000K if we use the collisional 
ionization curves of Benjamin & Shapiro (1996) for a gas that is cooling isobarically. (A similar 
argument arises from an upper limit for N V/N II, but the resulting constraint on the temperature 
is weaker.) There is considerably more Si III than Si II in the high velocity gas (Cowie, Songaila, 
& York 1979), but this is may be due to photoionization. Thus, to derive a lower limit for 
the temperature of the gas, we must use a typical equilibrium temperature for an H II region, 
somewhere in the range 8, 000 < T < 12, 000K (Osterbrock 1989). 



5.2.2. Immediate Postshock Gas 

Figure || shows that there is a broad absorption from C IV centered at a velocity of about 
—36 km s _1 , in addition to a narrower peak at about +20 km s -1 . The equivalent widths of 50 
and 22 mA for the broad, negative velocity components for the transitions at 1548.2 and 1550. 8A, 
respectively, indicate that logA r (C IV) = 13.1, a value that is in conflict with the upper limit 
logA^(C IV) = 11.9 obtained by Cowie, et al. (1979). Absorptions by Si III (1206.5A) and Al III 



Exposure identification: Z165040DM. 

Again, a medium resolution GHRS spectrum: Exposure identification: Z1650307T. 
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Fig. 8. — A medium resolution (R = 20,000) recording of the C IV doublet in a spectrum of £ Ori A 
taken with the G160M grating of GHRS on the Hubble Space Telescope. A correction of +0.06A 
has been added to the wavelength scale of calibrated file for exposure Z1650307T to reflect the 
offset in a 1554.9285A calibration line that may be seen in exposure Z1650306T. 
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(1854. 7 A) are also evident at —20 and —15 km s _1 , respectively] 1 " 1 ] We propose that these high 
ionization components arise from collisionally ionized gas behind the shock front. (Ultraviolet 
radiation from the shock front also helps to increase the ionization of the downstream gas.) The 
width of the C IV feature shown in Fig. [8| reflects the effects of thermal doppler broadening, 
instrumental smearing, and the change in velocity as the gas cools to the lowest temperature that 
holds any appreciable C IV. 



5.2.3. Properties of the Shock 

We return to our conjecture that the preshock gas flow is being intercepted by an obstacle at 
v « km s _1 , and thus the front itself is at this velocity. (While this assumption is not backed 
up by independent evidence, it is nevertheless a basic premise behind our relating the atomic 



absorption line data to our interpretation in §5.3 and §5.4 of how the H2 in Component 1 is formed 



in a region where there is a large compression and a temperature that is considerably lower than 
that of the immediate postshock gas.) The fact that the C IV feature does not appear at 74 times 
that of the high velocity (preshock) C II and N II features indicates that the compression ratio 
is less than the value 4.0 for strong shocks with an adiabatic index 7 = 5/3. This is probably a 
consequence of either the ordinary or Alfven Mach numbers (or both) not being very high. For 
example, if the preshock magnetic field and density were 5//G and uq = 0.1 cm -3 , the Alfven 
speed would be 29 km s _1 . For T = 20, 000K, the ordinary sound speed would be 21 km s _1 , and 
under these conditions the compression ratio would be only 2.67 [cf. Eq. 2.19 of Draine & McKee 
(1993)] if the magnetic field lines are perpendicular to the shock normal. This value is close to the 
ratio of velocities of the preshock and postshock components, (—94 km s _1 )/(— 36 km s _1 ) = 2.6. 
The immediate postshock temperature would be about 2.3 x 10 5 K. 

Our simple picture of a shock that is moderated by a transverse magnetic field adequately 
explains the velocity difference between the two atomic components, but it fails when we try to fit 
the kinematics of the much cooler gas where we find H2 . If we follow the material in the postshock 
flow to the point that radiative cooling has lowered the temperature to that of the preshock 
gas or below, we expect to have a final compression ratio equal to 3.7, i.e., the number that we 
would expect for an "isothermal shock" [cf. Eq. 2.27 of Draine & McKee (1993)]. This limited 
amount of compression would mean that the cool, H2-bearing gas would appear at a velocity of 
(—94 km s _1 )/3.7 = —25 km s _1 , a value that is clearly inconsistent with what we observe. 

A resolution of the inconsistency between the kinematics noted above and the theoretical 
picture of a shock dominated by magnetic pressure could be obtained if, instead of having the 
initial magnetic field lines perpendicular to the shock normal, the field orientation is nearly parallel 
to the direction of the flow. (Intuitively, this arrangement seems more plausible, since the field 



"Archive exposure identifications: Z165040CM and Z165040DM. 
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lines are likely to be dragged along by the gas.) The picture than can then evolve to the more 
complex situation where there is a "switch-on" shock, giving an initial moderate compression and 
a sudden deflection of the velocity flow and direction of the field lines. As described by Spitzer 
(1990a, b), this phase may then be followed by a downstream "switch-off" shock that redirects the 
flow and field lines to be perpendicular to the front and allows further compression of the gas up 
to values equal to the square of the shock's ordinary Mach number, i.e., the compression produced 
by a strong shock without a magnetic field. 

In order to obtain a solution for a switch-on shock, one must satisfy the constraint that 
the Alfven speed must be greater than slightly more than half of the shock speed [cf. Eq. 2.21 
of Draine & McKee (1993)]. Thus, we must at least double the Alfven speed of the previous 
example by either raising the preshock magnetic field, lowering the density, or both. If this speed 
equalled 58 km s _1 , the compression ratio in the switch-on region should be (94/58) 2 = 2.6, i.e., 
the square of the Alfven Mach number, a value that is again very close to our observed ratio of gas 
velocities on either side of the front. [There is a complication in deriving a compression ratio from 
an observation taken at some arbitrary viewing direction through a switch-on shock. Behind the 
front, the gas acquires a velocity vector component that is parallel to the front. For an inclined 
line of sight, this component can either add to or subtract from the projection of the component 
perpendicular to the front, which is the quantity that must be compared to the (again projected) 
preshock velocity vector when one wants to obtain a compression ratio. However in our situation 
it seems reasonable to suppose that a wind from £ Ori is the ultimate source of high velocity 
gas, and this in turn implies that the shock front is likely to be nearly perpendicular to the line 
of sight.] While this picture is still rather speculative, we will adopt the view that, through the 
mechanism of the switch-off shock, the magnetic fields do not play a significant role in limiting the 
amount of compression at the low temperatures where H2 could form. 

One additional piece of information is a limit on the preshock density uq. Cowie, et al. 
(1979) obtained an upper limit for the electron density n c < 0.3 cm' 3 from the lack of a detectable 
absorption feature from C II in an excited fine-structure level (assuming T = 10 4 K). Since there is 
virtually no absorption seen for lines of N I or O I at the high velocities in front of the shock, we 
can be confident that the hydrogen is almost fully ionized and thus the limit for n e applies to the 
total density. For the purposes of argument in the discussions that follow, we shall adopt a value 
no = 0.1 cm' 3 , as we have done earlier. 

5.3. Formation of H2 in a Warm Zone 

5.3.1. Reactions and their Rate Constants 

In the light of evidence from the atomic lines that a standing shock may be present, we move 
on to explore in a semiquantitative way the prospects that H2 forming behind this front could 
explain our observations. For several reasons, we expect that an initial zone where T £ 10 4 K will 
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produce no appreciable H2. At these temperatures the gas is mostly ionized, and for T > 18,000K 
collisions with electrons will dissociate H2 very rapidly (Draine & Bertoldi, in preparation). 
Furthermore, the column density of material at T > 6500K is not large because the cooling rate is 
high. As soon as the gas has reached 6500K, there is a significant, abrupt reduction in the cooling 
rate while there is still some heating of the gas by ionizing radiation produced by the much hotter, 
upstream material. These effects create a plateau in the general decrease of temperature with 
postshock distance [see Fig. 3 of Shull & McKee (1979)]. 

The 6500K plateau, extending over a length of approximately 2 x 10 16 nQ 1 cm, seems 
to be a favorable location for synthesizing the initial contribution of H2 that we could be 
viewing in the upper J levels. Its velocity with respect to much cooler gas should be about 
(—36 km s _1 ) x (6500K)/(2.3 x 10 5 K) = —1.0 km s _1 if the conditions are approximately isobaric. 
This velocity difference is consistent, to within the observational errors, with the shift between 
the peaks at J = 5 and J = 0, with the latter emphasizing molecules in the material that has 
cooled much further and come nearly to a halt. Considering that the fractional ionization over the 
temperature range 6500 > T > 2000K is 0.5 £ n e /n H £ 0.03 (Shull & McKee 1979), we anticipate 
that potentially important sources of H2 arise from either the formation of a negative hydrogen 
ion, 

H + e -> H~ + hv 

C 6 = 1.0 x 10~ 15 r 3 exp(-r 3 /7) cmV 1 (6) 

(T3 is the gas's temperature in units of 10 3 K) followed by the associative detachment, 

H~+H -► H 2 + e 

C 7 = 1.3 x lO^cmV 1 (7) 

or the production of H^" by radiative association, 

H + H+ -> H+ + hv 

C 8 = 4.1 x lO^cmV 1 (8) 

followed by its reaction with neutral atoms, 

H+ + H -> H 2 + H+ 

C 9 = 1.0 x 10~ 10 cm 3 s _1 (9) 

(Black 1978; Black, Porter, & Dalgarno 1981). The rate constants for the above reactions (plus 
the destruction reactions [l^ and [l3| below) are the same as those adopted by Culhane & McCray 
(1995) in their study of H2 production in a supernova envelope. Later, as the gas becomes cooler, 
denser and mostly neutral, we expect that the formation of H2 on the surfaces of dust grains, 

2H + grain — ► H2 + grain 

io _16 r - 5 

[applicable to n(H) 2 ] C 10 = - - - ^ 5 - ^ - 2 cmV 1 (10) 
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should start to become more important (Hollenbach & McKee 1979). We will address this 
possibility in § |5.4j . 

In order to evaluate the effectiveness of reactions ^ and ^ in producing H2 in the warm gas, 
we must consider the most important destruction processes that counteract the production of 
the feedstocks H~ (reaction ^) and (reaction ||). Radiative dissociation of H~ by uv starlight 
photons (i.e., the reverse of reaction ||), 

H + hv -> H + e 

fa = 1.9 x KnV 1 (11) 

is generally the most important mechanism for limiting the eventual production of H2 in partially 
ionized regions of the interstellar medium. The value for /^H] is adopted from an estimate for this 
rate of destruction in our part of the Galaxy by Fitzpatrick & Spitzer (1994). Less important 
ways of destroying H~ include recombination with protons, 

H~ + H+ -► 2H 

C12 = 7 x KT^^cmV 1 (12) 

and, of course, the production of H2 (reaction |7|). is destroyed by the reaction with electrons, 

H+ + e -> 2H 

C13 = 1.4 x 10- 7 r 3 -°- 4 cin 3 s- 1 (13) 

and the creation of H2 in reaction |9[ We can safely disregard the interaction of with H2 , 

+ H2 — ► + H 
C M = 2.1 x lO^cmV 1 (14) 

because C|fJn(H2) <C CJTJn(e) + CQn(H). Finally, our end product H2 is destroyed by 
photodissociation, 

U 2 + hv -> 2H 

(optically thin) (3 15 = 3.4 x 10 -11 s _1 . (15) 

Our adopted general value for makes use of Jura's (1974) calculation of /5|TJ = 5.4 x 10 _11 s _1 
for a flux of ^ttJ\ = 2.4 x 10 _6 erg cm _2 s _1 A 1 at 1000A, but rescaled to a local flux of 
1.5 x 10~ 6 erg cm _2 s _1 A calculated by Mezger, et al. (1982). A large fraction of this 
background may come from sources that are behind or within cloud complexes containing H2. 
If this is true, the stellar radiation in the cores of the most important Werner and Lyman lines 
is converted to radiation at other wavelengths via fluorescence (Black & van Dishoeck 1987), 
leading to a lower value for /3JTJ. The reduction in /3jTJ caused by self shielding of material within 
Component 1 is small: for log iV(H2) = 14.52 and 6 = 3 km s _1 it is only 32% (Draine & Bertoldi 
1996). We will discuss in § |5.3.3| how much the photodissociation of H2 could be increased by the 
gas's proximity to the hot, bright stars in the Orion association. 
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5.3.2. Expected Amount of H2 

We now investigate whether or not it is plausible that the above reactions can produce the 
approximate order of magnitude of H2 that we observe in the higher J levels of Component 1. For 
the condition that the preshock density no = 0.1 cm" 1 (§|5]1), we expect that the time scale for 
perceptible changes in temperature and ionization when 6500 > T > 2000K is about 4 x 10 11 s, 
a value that is much greater than the equilibration time scales j3^ for the production of H2, 
[t^+ C|7Jn(H + ) + Cln(H)]- 1 for H~, or [Cgn(e) + ^(H)]" 1 for H^. Thus, for a total density 
Tin = n(H + ) + n(H) and a fractional ionization / = n(H + )/nn the density of H2 at any particular 
location is given by a straightforward equilibrium equation 

n(H 2 ) = /(l - f) 2 n 2 u [F(R-) + F(H+) + F (grain) ]//? 15 (16a) 
F(R-) = ^31 (16b) 

F ^ = W^¥W) ' (16c) 



and 



F(grain) = C w /f (16d) 



In order to make an initial estimate for the amount of H2 that could arise from the warm, 
partly ionized gas, we must evaluate the integral of the right-hand side of Eq. fL6a| through the 
relevant part of the cooling, postshock flow. The structure of this region is dependent on several 



parameters that are poorly known and whose effects will be discussed in § 5.3.3 . As a starting 
point, however, we can define a template for the behavior of /, njj and T with distance by adopting 
the information displayed by Shull & McKee (1979) for a 100 km s _1 shock with no = 10 cnr 3 
and solar abundances for the heavy elements (their Model E displayed in Fig. 3). To convert to 
our assumed no = 0.1 cm -3 , we scale their densities n(H) and n(H + ) down by a factor of 100 and 
the distance scale up by the same factor. 

Over all temperatures, we discover that F^H^) is at least 100 times smaller than F(H~), 
and hence this term is not significant for our result, ingrain) is negligible compared to .F(H - ) 
at high temperatures, but its importance increases as the temperatures decrease: the two terms 
equal each other at T = 2000K, and F(grain) = 3.7F(H") at 1000K. Within the F(B~) term, 
the terms for photodestruction and recombination with H + in the denominator are about equal at 
T = 6500K, but the photodestruction becomes much more important at lower temperatures. 



The integral of the predicted n(H2) (Eq. |l6a| ) over a path that extends down to T = 2000K 
equals 3.6 x 10 13 cm~ 2 . This value is substantially lower than the amount of H2 that we observed 
in the higher J levels in Component 1 Ej= 2 N(H.2) = 2 x 10 14 cm -2 ; see Table [|. 
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5.3.3. Ways to Reconcile the Expected and Observed H2 

There are several effects that can cause significant deviations from the simple prediction 
for iV(H2) given above. First, if we accept the principle that the origin of the preshock flow at 
—94 km s _1 is from either a stellar wind produced by £ Ori (plus perhaps other stars in the 
association) or some explosive event in Orion, we must then acknowledge that the H2 production 
zone is probably not very distant from this group of stars that produce a very strong uv flux. 
As a consequence, we must anticipate that /5|7J could be increased far above that for the general 



interstellar medium given in Eq. 15. Eq. |16a| shows that this will give a reduction in the expected 
yield of H2 in direct proportion to this increase. (For a given enhancement of /37J, we expect that 
the increase in /5|7J will be very much less because the cross section for this process is primarily 
in the visible part of the spectrum where the contrast above the general background is relatively 
small.) Working in the opposite direction, however, is the fact that the stars' Lyman limit fluxes 
will supplement the ionizing radiation produced by the hot part of the shock front, thus providing 
heating and photoionization rates above those given in the model. The resulting higher level of 
/ and the increase of the length of the warm gas zone will result in an increase in the expected 
N(R 2 ). 

To see how important these effects might be, we can make some crude estimates for the 
relevant increases in the uv fluxes. In the vicinity of 1000A i.e., the spectral region containing 
the most important transitions that ultimately result in photodissociation of H2, the fluxes 
from e and a Ori at the Earth are 5.6 x 1CP 8 and 1.8 x 10 _8 erg cm _2 s _1 A , respectively 
(Holberg et al. 1982). We can assume that other very luminous stars that might make important 
contributions, such as 5, k and 1 Ori, have uv fluxes consistent with that of e Ori after a 
scaling according to the differences in visual magnitudes. The probable distance of the H2 from 
the stars is probably somewhere in the range 60 to 140 pc, as indicated by various measures 
of the transverse dimensions of shell-like structures seen around the Orion association (Goudis 
1982) (and assuming that the Orion association is at a distance of 450 pc from us). If we 
compare the far-uv extinction differences for 5 and e Ori reported by Jenkins, Savage & Spitzer 
(1986) to these stars' color excesses E(B— V) = 0.075, we infer from the uv extinction formulae 
of Cardelli, Clayton & Mathis (1989) that i?v = 4-6 and, again using their formulae, that 
A ° = 0.96 Mae. In the absence of such extinction, these plus the other stars should produce a 

1000A fa 1 

net flux ^ looa £ = 1-0 x 10~ 5 r{^Q erg cm _2 s _1 A , where rioo is the distance away from the stars 
divided by 100 pc. With rioo = 1) is enhanced over the value in Eq. [l5| by a factor of 7. 

Stars in the Orion association produce about 3.8 x 10 49 Lyman limit photons s _1 , and only a 
small fraction of this flux is consumed by the ionization of hydrogen in the immediate vicinity of 
the stars (Reynolds & Ogden 1979). From this estimate, one may conclude that the ionizing flux 
of ~ 10 6 photons cm^V 1 radiated by the immediate postshock gas (Shull & McKee 1979) could 
be enhanced by a factor approaching 30r{^Q, thus increasing the thickness of the region over which 
there is a significant degree of ionization and heating. 
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Another parameter that can influence the length of the zone where reactions ^ and || are 
important is the relative abundances of heavy elements. Here, the cooling is almost entirely from 
the radiation of energy by forbidden, semi-forbidden and fine-structure lines from metals - see 
Fig. 2 of Shull & McKee (1979). If these elements are depleted below the solar abundance ratio 
because of grain formation, the length of the warm H2 production zone must increase (Shull & 
Draine 1987). It is unlikely that the grains will been completely destroyed as they passed through 
a 90 km s _1 shock (Jones et al. 1994). 

Finally, it is important to realize that the outcome for N(H.2) should scale roughly in 
proportion to n$. The reason for this is that over most of the path, we found that /^j/nn was the 
most important term within denominator of the dominant production factor F(H~). This in turn 
makes n(H-2) scale in proportion to almost everywhere (note that / n^dl does not vary with 
no). 



5.3.4- Independent Information from an Observation of Si IF 

It is important to look for other absorption line data that can help to narrow the uncertainties 
in the key parameters discussed above. One such indicator is the column density of ionized silicon 
in an excited fine-structure level of its ground electronic state (denoted as Si II*). This excited 
level is populated by collisions with electrons, and the balance of this excitation with the level's 
radiative decay (and collisional de-excitations) results in a fractional abundance 

log ( ^n? ) = log "(e) " 0.5 log T 3 - 2.54 (17) 

(Keenan et al. 1985). In conditions where the hydrogen is only partially ionized, we expect that 
ra(Si ++ )/ra(Si + ) will be much less than ra(H + )/re(H) because ionized Si has a larger recombination 
coefficient (Aldrovandi & Pequignot 1973) and a smaller photoionization cross section (Reilman 
& Manson 1979) (its ionization potential of 16.34 eV is also greater than that of hydrogen). Thus, 
for situations where / is not very near 1.0, it is reasonably safe to assume that virtually all of the 
Si is singly ionized. If, for the moment, we also assume that the Si to H abundance ratio is equal 
to the solar value, we expect that 

n(Si II*) = 10- 7 /n|T 3 ~ 1/2 cm" 3 . (18) 

As we did for H2, we can integrate the expression for n(Si II*) through the modeled cooling zone 
to find an expectation for the column density iV(Si II*) = 5.8 x 10 em~ 2 . 

A very weak absorption feature caused by Si II* at approximately the same velocity as our 
Component 1 can be seen in a medium resolution HST spectrum^ of £ Ori A that covers the very 
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strong transition at 1264. 730A. Our measurement of this line's equivalent width was 13.4 ± 3mA, 
leading to iV(Si II*) = 1.0 x 10 12 cm -2 , a resuhf^| that is almost twice the prediction stated above. 

From our result for Si II*, we conclude that the combined effect of the stars' ionizing flux 
and a possible increase in no over our assumed value of 0.1 cm -3 could raise / n(e)nndl by not 
much more than a factor of two. However, we have no sensitivity to the possibility that metals are 
depleted since the decrease in the abundance of Si would be approximately compensated by the 
increase in the characteristic length for the zone to cool (assuming the primary coolants and Si are 
depleted by about the same amount). Thus, it is still possible that the our calculation based on a 
model with solar abundances will result in an inappropriate (i.e., too low) value for the expected 
iV(H 2 ). 

5.3.5. Coupling of the Rotational Temperature to Collisions 

One remaining task is to establish that the conditions in the H2-formation zone are such that 
collisional excitation of the higher J levels can overcome the tendency for the molecules to move 
to other states through either radiative decay or the absorption of uv photons. Tawara et al. 
(1990) summarize the collision cross sections as a function of energy for excitations J = —* 2 and 
J = 1 — ► 3 by electrons. We calculate that these cross sections should give a rate constant of about 
1 x 10~ 10 T 3 3/2 cm s for T3 ;> 2. Thus, for n e ^ 1 cm and T3 rj 5 the collisions can dominate 
over radiative transition rates of about 3 x 10 -10 s _1 for J = 2 and 3. To collisionally populate 
J = 5 which can decay at a rate of 1 x 10 -8 s _1 to J = 3, we would need to have n e ^ 10 cnr 3 
just to match the radiative rate, assuming that the collisional rate constant is not significantly 
lower than what we calculated for J = 2 and 3. 

5.4. Further Formation of H2 in a Cool Zone 

Additional formation of H2 molecules probably takes place in gas that has cooled well below 
2000K and is nearly fully recombined. We are unable to distinguish between this gas and the 
material that was originally present as an obstruction to the high velocity flow to create the 
bow shock. To obtain an approximate measure of the total amount of cool, mostly neutral gas, 
we determined / N a {y)dv over the velocity interval — 10 < v < +5 km s _1 for the N I line at 
1134. 165A which does not appear to be very strongly saturated. The relative ionization of nitrogen 
should be close to that of hydrogen (Butler & Dalgarno 1979), and this element is not strongly 
depleted in the interstellar medium (Hibbert, Dufton, & Keenan 1985). Our conclusion that 



13 From the equivalent width of 3.4mA (no error stated) for the 1194. 49A line of Si II* reported by Drake & Pottasch 
(1977), one obtains a somewhat lower value, logiV(Si IF) = 11.6 
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logiV(N I) = 14.78 leads to an inferred value for log iV(H) equal to 18.73. According to a modelf^l 
for a 90 km s" 1 shock of Shull & McKee (1979), logiV(H) = 18.40 is the amount of material that 
accumulates by the time T reaches 1000K. Hence, from our measure of the total N(H) (but indeed 
an approximate one), we estimate that the amount of gas at T < 1000K is about comparable to 
that at the higher temperatures. 

An insight on the conditions in the cool, neutral zone is provided by the populations 
of excited fine-structure levels of C I. Jenkins & Shaya (1979) found that logp/k = 4.1 
in the part of our Component 1 that carries most of the neutral carbon atoms. If we take 
as a representative temperature T = 300K, the local density should be nn ~ 40 cm~ 3 
and C|io| = 1.8 x 10 cm 3 s -1 . With the general interstellar value for /5JTJ, we expect an 
equilibrium concentration n(H.2)/nn = 2 x 10~ 5 . When we multiply this number by our estimate 
iV(H) = 3 x 10 18 cm -2 , we find that we should expect to observe logiV(H2) = 13.8, a value that, 
after considering the crudeness of our calculations, is acceptably close to our actual measurements 
of H2 in J = and 1, the levels that arise primarily from the coolest gas. If /3^| is significantly 
enhanced by radiation from the Orion stars (§ 5.3.3]) , we would then have difficulty explaining the 
observations. 



6. Summary 

We have observed over 50 absorption features in the Lyman and Werner bands of H2 in the 
uv spectrum of £ Ori A. An important aspect of our spectrum is that it had sufficient resolution 
to detect in one of the velocity components (our Component 1) some important changes in the 
one-dimensional velocity distributions of the molecules with changing rotational excitation J. 
The main focus of our investigation has been to find an explanation for this result, since it is a 
departure from the usual expectation that the rotational excitation comes from uv pumping, an 
effect that would make the profiles look identical. A smaller amount of broadening for higher J is 
also seen for Component 3, a component that has much more H2 than Component 1. 

In Component 1, we have found that as J increases from to 5 there is a steady increase 
in the width of the velocity profile, combined with a small drift of the profile's center toward 
more negative velocities. We have shown that the pumping lines are not strong enough to make 
differential shielding in the line cores a satisfactory explanation for the apparent broadening of 
the J levels that are populated by such pumping. While one might resort to an explanation that 
unseen, additional H2 could be shielding light from a uv source (or sources) other than £ Ori, 
we feel that this interpretation is implausible, and, moreover, it does not adequately explain the 
differences that we see between the profiles of J = and 1. 



1 This model is not exactly applicable to our situation, since it has a compression ratio of 4 instead of our value of 
2.6 
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One could always argue that the absorption that we identify as Component 1 is really a 
chance superposition of two, physically unrelated regions that have different rotation temperatures 
and central velocities]^] While not impossible, this interpretation is unattractive. It requires a 
nearly exact coincidence of the two regions' velocity centers to make up a component that stands 
out from the rest of the H2 absorption and, at the same time, shows smoothly changing properties 
with J at our velocity resolution. We feel that the most acceptable interpretation is the existence 
of a coherent region of gas that, for some particular reason that has a rational explanation, has 
systematic changes in the properties of the material within it that could produce the effects that 
we observe. One phenomenon that fits this picture is the organized change in temperature and 
velocity for gas that is cooling and recombining in the flow behind a shock front. The excess width 
of the higher J lines could arise from both higher kinetic temperatures and some velocity shear 
caused by the steady compression of the gas as it cools. 

Our concept of a shock is supported by evidence from atomic absorption lines in the spectrum 
of C Ori A. We see features that we can identify with both the preshock medium and the immediate 
postshock gas that is very hot. If this interpretation is correct (and not a misguided attempt to 
assign a significant relationship between atomic components with different levels of ionization at 
very different velocities), we can use the atomic features to learn much more about the shock's 
general properties. 

We start with the expectation that the coolest molecular material, that which shows up in the 
lowest J levels at v = —1 km s _1 , is in a region containing gas that is very strongly compressed 
and thus nearly at rest with respect to the shock front. The atomic features of C II and N II that 
we identify with the preshock flow appear at a velocity of —94 km s _1 with respect to the cool H2. 
Hence this is the value that we normally associate with the "shock velocity." This preshock gas 
also shows up in the lines of C III, N III, Si III, S III and Al III. An upper limit to its temperature 
T < 20, 000K results from the apparent lack of C IV that would arise from collisional ionization at 
slightly higher temperatures. The temperature could be as low as typical H II region temperatures 
(8, 000 < T < 12, 000K) if uv photons are the main source of ionization. 

Absorption features from more highly ionized gas at around —36 km s _1 that show up in the 
C IV doublet indicate that the initial compression factor is only 2.6, a value that is significantly 
lower than the usual 4.0 expected for a shock with a high Mach number. Reasonable numbers for 
the preshock density, temperature and magnetic field strength (0.1 cm -3 , 20,000K and 10/uG) can 
explain this lower compression factor and establish a switch-on shock. However, except for a limit 
no < 0.3 cm~ 3 that comes from the lack of C II* absorption, we have no independent information 
that can distinguish between these somewhat arbitrary assignments and other, equally acceptable 



15 There is a good way to illustrate how this creates the effect that we see in Component f. Imagine that we 
recorded the H2 lines at a resolution that was so low that Complexes I and 3 were not quite resolved from each other. 
We would then see features that got broader with increasing J, and their centers would shift toward the left. This is 
qualitatively exactly the same effect that we see on a much smaller velocity scale within Component 1. 
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combinations. 

To overcome the problem that there seems to be a low initial compression of the gas but 
eventually the densities are allowed to increase to the point that H2 can form, we invoke the 
concept of an oblique magnetic shock, where the theoretical models outline the existence of two 
discontinuities, a "switch on" front and a "switch off" front. However, we do not attempt to 
explore the validity of this picture in any detail. 

Neglecting complications that are introduced by the oblique shock picture, we expect that as 
the gas flow cools to temperatures significantly below the immediate postshock temperature, it 
decelerates and begins to show ions that have ionization potentials below that of C IV, such as 
Si III (at —20 km s _1 ) and Al III (—15 km s _1 ). At temperatures somewhat below 10 4 K, the gas 
should be still partially ionized and at a density bh > 60no. These conditions favor the production 
of H2 via the formation of H~ and its subsequent reaction with H to produce the molecule plus an 
electron, rather than the usual formation on grains that dominates in cool clouds. Our observation 
of the Si II* absorption feature at a velocity near km s _1 indicates that it is unlikely that 
no -C 0.1 cm~ 3 , and thus the density in the molecule forming region is high enough to insure that 
the photodetachment of H~ does not deplete this feed material to the point that the expected 
abundance of H2 is well below the amount that we observe. 

It is possible that the uv flux from the Orion stars could enhance the H2 photodissociation 
rate /^j to a level that is far above that which applies to the average level in our part of the Galaxy. 
If this is true, then we must make a downward revision to our prediction that log N(H.2) = 13.56. 
At the same time, however, additional ionizing photons from the stars could lengthen the warm 
H2 production zone, and this effect may gain back a large amount of the lost H2. A reduction 
in the metal abundance in the gas may also lengthen the zone, giving a further increase in the 
expected H2. 

As the gas compresses further and becomes almost fully neutral, the H~ production must 
yield to grain surface reactions as the most important source of molecules. Using information from 
the C I fine-structure excitation, we can infer that the density in the cool gas is sufficient to give 
n(H2)/nH equal to about half of what we observed, if we assume that most of the H2 in the J = 
and 1 states comes from the cool region. 

On the basis of a diverse collection of evidence and some rough quantitative calculations, 
we have synthesized a general description of the cooling gas behind the shock and have shown 
that H2 production within it could plausibly explain the unusual behavior in the profiles that 
we observed. Obviously, if one had the benefit of detailed shock models that incorporated the 
relevant magnetohydrodynamic, atomic and molecular physics, it would be possible to substantiate 
this picture (or perhaps uncover some inconsistencies?) and narrow the uncertainties in various 
key parameters. Also, more detailed models should allow one to address certain questions that 
are more difficult to answer, such as whether or not more complex chemical reactions play an 
important role in modifying the production of H2; we have identified only a few good prospects. 
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For instance, is there enough Ly-a radiation produced in the front (or in the H II region ahead of 
it) to make the formation by excited atom radiative association (i.e., H(n = 1) + H(n = 2) — ► H2 
+ hf) an important additional production route (Latter & Black 1991)? On the observational 
side, we expect to see very soon a vast improvement in the amount and quality of data on atomic 
absorption lines toward ( Ori A. Very recently, the GHRS echelle spectrograph on HST obtained 
observations of various atomic lines at extraordinarily good resolution and S/N. 

Support for flying IMAPS on the ORFEUS-SPAS-I mission and the research reported here 
came from NASA Grant NAG5-616 to Princeton University. The ORFEUS-SPAS project was 
a joint undertaking of the US and German space agencies, NASA and DARA. The successful 
execution of our observations was the product of efforts over many years by engineering teams at 
Princeton University Observatory, Ball Aerospace Systems Group (the industrial subcontractor 
for the IMAPS instrument) and Daimler-Benz Aerospace (the German firm that built the 
ASTRO-SPAS spacecraft and conducted mission operations). Most of the development of the 
data reduction software was done by EBJ shortly after the mission, while he was supported by 
a research award for senior U.S. scientists from the Alexander von Humboldt Foundation and 
was a guest at the Institut fur Astronomie und Astrophysik in Tubingen. We are grateful to 
B. T. Draine for valuable advice about the different alternatives for interstellar shocks. B. T. 
Draine, L. Spitzer, and J. H. Black supplied useful comments on an early draft of this paper. 
Some of the conclusions about atomic absorption features are based on observations made with 
the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by AURA under NASA contract NAS 5-26555. 
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